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Abstract 

A superconductive model characterized by a third order parabolic operator 
C £ is analysed. When the viscous terms, represented by higher - order deriva- 
tives, tend to zero, a hyperbolic operator Co appears. Furthermore, if V £ is the 
Dirichlet initial boundary - value problem for C e , when C e turns into Co, V e 
turns into a problem Vo with the same initial - boundary conditions as V £ . The 
^ . solution of the nonlinear problem related to the remainder term r is achieved, 

as long as the higher-order derivatives of the solution of Vo are bounded. More- 
c/2 . over, some classes of explicit solutions related to Vo arc determined, proving the 

existence of at least one motion whose derivatives are bounded. The estimate 
shows that the diffusion effects are bounded even when time tends to infinity. 
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1 Introduction 

We deal with the following third order parabolic equation: 

(d xx - Xd x ) (eut + u) - d t {u t + au) = F(x,t,u) (1) 

which characterizes a lot of phenomena (see, for instance, [1] and references therein). 
In particular, when F = sinu — 7 (7 = const), the semilinear equation (1) describes 
the evolution of the phase u inside an exponentially shaped Josephson junction 
(ESJJ) [2,3]. 

The (ESJJ) provides several advantages in comparison with rectangular junc- 
tions and the related bibliography is large and extensive. [1-3, 16-21]. Besides, 
interesting applications of Josephson junctions are SQUIDs which can be applied in 
a lot of fields [4], Indeed, they are very versatile and they are used in medicine for 
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the Magnetoencephalography (MEG) and for Magnetic resonance imaging (MRI) or 
in cosmology, beeing SQUIDs extremely sensitive detectors of magnetic and electric 
signals from any origin. Moreover, in geophysics they are used as magnetometers 
to discover mineral deposits, or as X-ray detectors to screen radioactive materi- 
als. [5,6]. Furthermore, they are suitable in the development of smart materials 
such as superconductor cables, used to improve power grids. [7,8] 
Equation (1) is equivalent to the equation: 

eA 2 eA 2 a a 

eu xxt ~ u tt + u xx - {a + — —)u t = —u+e~2 x F(x,t,e2 x u) (2) 



and from this, letting 



4 



eA 2 1 , , A 2 e a 

a + and 



4 e 4 e 

it is possible to obtain the following system : 



du d 2 



v —au — ee a x F(x, t, e 2 x u) 



dt dx 2 

d v 1 a 

— = bu v - e u t F t {x,t,e2 x u) 

at e 



(3) 



which is related to biological dynamical systems such as, for instance, the Fitzhugh 
Nagumo model (FHN) which reproduces the neurocomputational features of the 
neuron [9, 10]. 

So, as it has already been underlined in [11-15] , by means of equation (1) it is 
possible to treat both superconductive problems and neuroscientific ones. 

In this paper, the Dirichlet initial boundary value problem for equation (1) is 
considered and, according to the theory of the phenomena related to (ESJJ), as 
source term, it is assumed that F = smu — 7 (7 = const). 

Equation (1) is a semilinear hyperbolic equation, perturbed by viscous terms 
which are described by higher - order derivatives with small diffusion coefficients e. 
When e = 0, the parabolic equation (1) turns into a hyperbolic equation: 

(d xx - Xd x )uq - d t (d t + a)u = F(x,t,u ), (4) 

and the influence of the dissipative terms e (u xx t — Xu x t) on the wave behaviour of 
no can be estimated when the difference 

r(x,t,e) = u(x, t, e) — uq(x, t) (5) 

is evaluated. 

A similar problem has already been studied in [22], when the coefficient A in (1) is 
equal to zero and for finite intervals of time. Moreover, in [23,24], for A = and a = 



2 



0, an asymptotic approximation is established by means of the two characteristic 
times: slow time r = et and fast time 6 = t/e. 

In this paper, in order to evaluate the diffusion effects, a rigorous estimate of the 
reimander term r, defined in (5), is done. Indeed, by means of a Fourier series with 
properties of rapid convergence, the solution in both linear and non linear case of 
the problem related to r, is determined. Furthermore, some classes of solutions of 
the hyperbolic equation have been determined, proving that there exists almost a 
solution whose derivatives are bounded. Finally, as soon as the hyperbolic equation 
admits solutions with limited derivatives, an estimate for the remainder term is 
achieved by proving that the diffusion effects are both bounded when e goes to zero 
and when the time t tends to infinity. 

2 Statement of the problem 

If T, I are two positive constants and 

n T = {(x, t) : < x < I, 0<t<T}, 
let us consider the initial boundary value problem V e : 

(d X x - A d x ) (eu t + u) - d t {u t + au) = F(x, t, u), (x, t) e Qt, 
< u(x, 0) = ho(x), ut(x, 0) = h±(x), x€[0,l], (6) 
u(0,t) = ifo, u(l,t) = <pi, <t <T. 

where F, hi, ifi{i = 0,1) are arbitrary specified functions and e, a, and A are 
positive constants. 

If e =0, the parabolic equation (6)i turns into the hyperbolic equation (4) and 
the new problem Vo for uo(x,t) has the same initial boundary conditions of V £ . 

To evaluate the influence of the dissipative terms, letting 

F(x,t,r) = F(x,t,r + u ) - F(x,t,u ) - ed xt (d x - X)u , (7) 

the following problem V related to the remainder term r, defined in (5) has to be 
considered: 

(d xx - \d x ) (ed t +l)r-d t (d + a)r = F(x,t,r), {x,t)€Qr, 

r(x,0) = 0, r t (x,0) = 0, x G [0,1], (8) 

r(0,t)=0, r(l,t)=0, 0<t<T. 

As for the model of superconductivity, the problem (8) is characterized by the 
source force: 
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F = sin (r + ii ) -sinu - ed xt (d x - \)u . 



0) 



3 The Green function and its properties 

Let us introduce the linear third order operator: 

Ce = (d xx - Xd x )(ed t + 1) - d t {d t + a). (10) 

The Green function of C £ has already been determined in [1] by means of Fourier 
method. Let 



7n = ^j; K = (t£ + A 2 /4), g n = -(a + sb n ), oj n = V dl ~ K (H) 

and 

G n {t) = — e~ 9nt sinh(uj n t), (12) 
standard techniques allow to obtain the Green function as: 

G(x, t, f) = - e2 x y~] G n (t) sinrf n ^ sin^ n x. (13) 
' n=l 

The following theorem, proved in [1], ensures that this series is characterized by 
properties of rapid convergence and that it is exponentially vanishing as t tends to 
zero. 

Theorem 1. Whatever the constants a, e, A may be in !R + , denoting by a\ = 
a + e A 2 /4 and with 

7T 2 ax + e{7r/l) 2 . 
P\ = — 2~, <?A = , S = min(p x ,q\), (14) 

i/ie function G(x,£,t) defined in (13) and all its time derivatives, are continuous 
functions in f2y and it results: 

f)jr< 

\G(x,£,t)\ < Me~ 5t , |^-| < N,e~ 5t , j e N (15) 

where M,Nj are constants depending on a, A,e. 
Moreover, one has: 

Idi'HeGt + G)\ <A e~ st , (i = 0,1,2) (16) 
where A{ (i = 0,1,2) are constants depending on a, e, A. □ 
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The uniform convergence proved in theorem 1 allows to prove that [1]: 

C e G = (d xx - Xd x )(eG t + G) - d t (G t + aG) = 0. (17) 

4 The remainder term 

To make an estimate of the solution to the problem V related to the remainder term 
r, two cases should be distinguished according to the source term being known or 
depending by the unknown function. 

Let us assume F(x,r,t) = f(x,t). In this circumstance the remainder term r 
can be expressed by an explicit form. Indeed, standard computations (see f.i. [25]) 
ensure that the solution of the linear problem (8) is given by 

r f (x,t)=- f dr f G(x,Z,t-T)f(Z,T)d£. (18) 
Jo J 

Otherwise, if the source term F is non linear, an integro differential formulation of 
the problem can be obtained. 
Let us introduce the following 

Conditions A: Let us assume that function F( x,t,u) is defined and continuous 
on the set 

Dt = { (x, t, u) : (x, t) 6 , —oo < u < oo } (19) 

and that it is Lipschitz continuous in u i.e. there exists a constant K such that 
the inequality 

\F{x,t,U\) - F(x,t,u 2 )\ < K \ui-u 2 \ (20) 

holds for all (u±, u 2 ) and (x, t) £ f^. 

Under Assumption A, by means of the fixed point theorem, it is possible to 
prove [26-28] that the non linear problem V admits a solution given by 

r(x,t)= f dr f G(x,C,t-r)F(C,T,r(C,r))dC (21) 
Jo Jo 

5 Explicit solutions of the hyperbolic equation 

Let us consider the semilinear second order equation: 

u ,xx - u 0)tt - au 0)t -Xu 0tX = sin-uo - 7 (22) 

When A = 0, (22) represents the perturbed sine - Gordon equation (PSGE) 
and there is plenty of literature about its classes of solution, (e.g. [29-31] and 
bibliography therein.) 
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Now, let V De an arbitrary function and let us consider II(^) as: 

II(V0 = 2 arctan (23) 



such that 



sin II('0) 



1 



cos n(^) = — tanh(^). 



(24) 



cosh(?/>) ' 

By means of function (23) it is possible to find a class of solutions of equation (22). 
Indeed, it is possible to verify that the following function: 



uo=2II[/(£)] with £ 
is a solution of (22) provided that one has: 



x — t 
a — A 



(A — a) UQ t t = sin-uo — 7. 



(25) 



(26) 



Moreover, since (24) and beeing IT 



cosh / ' 



(A - a) u ,t 



sinuo = 2 sin II cos II 



it results: 



cosh / 



tanh / 



cosh / 

so, from (26), one deduces that function / must satisfy the following equation: 

df 



tanh / +7/2 cosh/ 



-de 



(27) 



Since physical experiences show that generally < 7 < 1, we point our attention to 
those cases in which it results: 



uq = 4arctan(y + \Jy 2 + 1 ). 
So, let h be an arbitrary constant of integration, one obtains: 



y 



y 



1 - (£ - h) 



l + {-h 
Moreover, assuming 7 < 1, let 



when 



when 



7 = 0, 



7 = 1. 



(28) 

(29) 
(30) 



a= ±v/r 



and 



(3 = 1 + A. 



For 7 2 / /3 2 , it results 
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1-he* 

Remark 1: In the case 7 = 0, if denoted by 



X — t 



z(x,t) = e~ ^ + \j e~ 2 ^ + 1, (32) 
equation (22) admits the solution 

uo(x,t) = 4 arctan [z(x, t)]. (33) 
and it is possible to prove that 

1 y(y 2 -i) 



UQ,xt(x,t) -- 



2{a-X) 2 (y 2 + l) 2 

1 y (y 6 - 5y 5 - 5y 2 + 1) 



(34) 



2(a-A) 3 (y 2 +1) 4 

which are bounded for all (x,t) € Qt- 

6 Estimates for the reimander term 

Let us assume e = 0, and let uq be a solution of the problem Vq. 

Theorem 2. If there exist two positive costards h, k depending only on a and A 
such that 

\uo,xt(x,t)\ < h; \u ,xxt{x,t)\ < k, (35) 

denoting by 



S{t) = sup \r(x,t,e)\, (36) 



the following inequality holds: 

0<S(t)< Ml{h s +k)£ e( 1 — (37) 
beeing M a positive constant depending on a, A,e and 5 is defined by (14)- 
Proof : Let us consider function F defined in (9): 

F = sin(r + u ) -sinu - ed xt (d x - A)n . (38) 

That source term satisfies conditions A and so, by means of hypotheses (35), it 
results: 
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\r(x,t,e)\ < f dr f G(x,£,t - r) [ |r(f,t,e) | + e(h + k)] d£. (39) 
Jo Jo 

By properties of the Green function G, and in particular since (15)i, one obtains: 

0<S(t)<Mlf e ~ 8{t ~ T) [5(r) + s(h + k)]dr (40) 
Jo 

So, by means of Gronwall Lemma, inequality (37) holds. □ 
Inequality (37) shows that, under assumptions of theorem 2, in the evolution of 
u(x,t) the correction due to the e - terms is bounded even when t tends to infinity, 
and it vanishes as soon as e tends to zero. 

Remark 2 Formulas (34) show that the class of functions satisfying hypoteses 
of Theorem 2 is not empty. 
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